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Abstract 

The quark Dyson-Schwinger equation and meson Bethe-Salpeter equation are stud- 
ied in a truncation scheme that extends the rainbow-ladder approximation such 
that, in the chiral hmit, the isovector, pseudoscalar meson remains massless. Quark- 
quark (diquark) correlations, which are bound in rainbow-ladder approximation, 
are destabilised by repulsive contributions that only appear at higher order in the 
Bethe-Salpeter kernel. The net effect of higher order terms on the meson bound- 
state masses is small. 
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1. Introduction. 

The analytic structure of quark and gluon propagators can provide an under- 
standing of confinement. Should these propagators have no Lehmann represen- 
tation then a sufficient condition for confinement is satisfied; i.e., the absence 
of quark and gluon production thresholds in iS-matrix elements describing 
colour-singlet to singlet transitions. This in itself, however, neither precludes 
nor entails the existence of bound-states, whether coloured or not. The exis- 
tence of two-body bound-states can only be established via the Bethe-Salpeter 
equation [BSE]. 

The homogeneous Bethe-Salpeter equation is derived under the assumption 
that the associated two-body T-matrix has a pole in a given channel. The ab- 
sence of a solution to the equation contradicts this assumption and establishes 
that no bound-state exists with the quantum numbers of the channel under 
consideration. 
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The intrinsically nonperturbative nature of bound-state problems, and the 
complex structure of the QCD vacuum, suggests that one should employ 
dressed gluon and quark propagators in constructing the kernel of the BSE. 
These can be obtained via Dyson-Schwinger equations [DSEs]. The Bethe- 
Salpeter equation is one type of DSE. 

There have been many studies of meson and diquark spectroscopy using this 
framework; a summary can be found in Ref. [1]. Typically, such studies em- 
ploy an Ansatz for the gluon propagator in solving the rainbow-approximation 
quark-DSE and pair the input gluon propagator and calculated quark propa- 
gator to construct the dressed-ladder approximation kernel for the meson/di- 
quark BSE. The BSE is then solved to obtain the spectrum. This rainbow- 
ladder truncation has the feature that Goldstone's theorem is manifest; i.e., 
in the chiral limit, when the current-quark mass niq — 0, the pion is a zero- 
mass bound-state of a strongly-dressed quark and antiquark [2]. With a few- 
parameter model for the gluon propagator, it can be used to provide a good 
description of the light-light, light-heavy and heavy-heavy meson spectrum 
and decays [3]. However, it has the defect that it admits (3)c-diquark bound- 
states [4]; such coloured states have not been observed. 

Herein we consider a truncation scheme that allows for a systematic improve- 
ment in the construction of the kernels of the quark-DSE and meson/diquark- 
BSE. This procedure ensures that the pion remains a Goldstone boson at 
every order. The first correction to ladder-rainbow approximation introduces 
a repulsive term in the Bethc-Salpeter kernel. Using a model gluon propaga- 
tor, singular in the infrared and with no Lehmann representation, we solve the 
quark DSE and obtain a quark propagator that also has no Lehmann represen- 
tation. Pairing these gluon and quark propagators in the meson channel, the 
repulsive term is almost cancelled by attractive terms of the same order and 
therefore the higher order terms only lead to a small change in meson masses. 
However, due to the algebra of SU{3)c, the repulsive term is stronger in the 
diquark channel and is not cancelled by the attractive terms. This entails that 
there is no stable diquark bound-state. 

2. Dyson-Schvi^inger and Bethe-Salpeter equations. 

In Euclidean metric, with 7^ hermitian and {7^,7;/} = 25^,^, the propagator 
for a dressed quark has the form S{p) = —ij ■ p<Jv{p'^) + crsip^)- It is obtained 
as a solution of the quark DSE 



S (p) = ij ■ p + niq + T,{p) (1) 
^ij.p + rriq + g^ I ^^^^ S{k) re(A;,p)y D^,{p - k), 

where niq is the current-quark mass, r^(/c,p) is the dressed quark-gluon vertex 
and D^i,{k) is the dressed gluon propagator, which, in Landau gauge, can be 
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written 



The mass and Bethe-Salpeter amplitude for a quark-antiquark bound-state, 
Fm, is obtained by solving 



^ffip^P)- (3) 
/ ^KZ'''''''{k,p;P) {S{k + \P)TM{k-P)S{k-\P)f'' , 

where KM{k,p; P) is the kernel, P is the total-momentum of the system, 
k,p are the internal and external relative quark-antiquark momenta and the 
superscripts are associated with the colour, flavour and Dirac structure of the 
amplitude; i.e., E = {icifyio}- 

The analogue of Eq. (3) for quark-quark systems is 



^d{p;P)- (4) 
/ {s(k+'^P)rUk;P)s^{-k+^,p)f" , 

where "T" denotes matrix-transpose. The absence of a solution to this equa- 
tion entails that diquarks do not appear in the strong-interaction spectrum. 

In the isovector channel, each contribution to K^^''^^{k,p; P) has a direct 
analogue in K^^''^^ {k,p; P); an analogue that can be obtained via the re- 
placement 

T 

(5) 

({A"/2}^^]^ are the generators of SU{3)c) in each antiquark segment of the 
meson kernel, which can be traced unambiguously from the external antiquark- 
line of the meson amplitude. This means that there are no contributions in 
the diquark channel that are not also present in the isovector meson channel. 

Rainbow-ladder approximation. 

The rainbow approximation to the quark DSE is defined by the choice r^(A;, p) = 
7^ in Eq. (1); the ladder-approximation to the meson BSE by 

P) [S{k + lP)rMik; P)S{k - \P)f" (6) 



l,^S{-k) 
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= -g^ D,,{p - k) (^7m Y S{k + iP)rM(fc; P)S{k - 1P)7. y 
in Eq. (3). 

In the chiral limit the rainbow-approximation DSE and the pseudoscalar- 
meson BSE in ladder-approximation are equivalent when = 0; i.e., one 
necessarily has a massless, pseudoscalar bound-state when chiral symmetry 
is dynamically broken[2]. Goldstone's theorem is manifest. In any truncation, 
such an identity between the quark DSE and the BSE in the isovector, pseu- 
doscalar, meson channel is sufficient to ensure that this meson is a Goldstone 
boson [5]. This provides for a straightforward understanding of the dichotomy 
of the pion as both a Goldstone boson and quark-antiquark bound-state. 

A gluon propagator that has no Lehmann representation may be interpreted 
as describing a confined particle. Such a gluon propagator can, via Eq. (1), 
yield a quark propagator with the same property [6]. Hence, in rainbow- ladder 
approximation, one can develop a spectroscopic model of mesons with confined 
quarks and gluons. 

However, at this order the diquark BSE, Eq. (4), admits solutions; i.e., one 
has (3)c diquark states in the spectrum [4], which are not observed. (It follows 
from the algebra of the generators of SU (3)c that, in this approximation, there 
are no (6)c diquark solutions just as there are no (8)c meson solutions.) 

3. Next-order Truncation 

As the next level of truncation for the quark DSE we consider Eq. (1) with 

n{k,p) = 7. + ^ / 9'D,,{P - I) Ip S{l + k- p)7. 5(07. ■ (7) 

This is the first-order correction of the vertex by the dressed-gluon propagator 
and the dressed-quark propagator, which is obtained as the solution of this 
DSE. Here we omit the exphcit 3— gluon vertex that could contribute at this 
order. This entails that we explore 3— and 4— body forces only to the ex- 
tent that they are incorporated via the nonperturbative dressing of the gluon 
propagator. 

For the purpose of illustration and clarity we employ a model gluon prop- 
agator [6] that allows for an algebraic solution of the Dyson-Schwinger and 
Bethe-Salpeter equations: 

A{e) = 16Tr^G5\k) , (8) 
in Eq. (2), where G — rf /A with 77 a mass scale. This one-parameter model 
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possesses the infrared enhancement, due to the 3— gluon vertex, suggested by 
the studies of Refs. [7] but underestimates the interaction strength for k"^ > 0. 
This is appropriate in a study whose focus is confinement and dynamical chiral 
symmetry breaking, which are infrared effects. The quahtative features of our 
results are not sensitive to this choice. 

Using Eqs. (7) and (8) in Eq. (1) one obtains^ 

Neglecting the O(G^) term yields the rainbow-approximation quark DSE. 
Goldstone Theorem. 

An important consideration in extending the kernel of the Bethe-Salpeter 
equation is to ensure that in doing so one preserves the Goldstone boson 
character of the pion. To this end we observe that the ladder- approximation 
kernel in the meson BSE can be obtained directly from the expression for 
in the rainbow-approximation quark DSE via the replacement 

7/. S{k) 7, -> 7^ S{k + P/2) TM{k, P) S{k - P/2) 7, , (10) 
which is illustrated in the top diagram of Fig. 1. 

In the isovector channel, this procedure can be implemented at every order; 
i.e., in every term of the quark DSE one may sequentially replace each explicit, 
internal quark propagator in this way, which is illustrated for Eq. (9) in Fig. 1. 
This generates all contributions of a given order to the kernel and ensures that 
Goldstone's theorem is preserved at that order, as will become clear below. 
Applying the procedure to vacuum polarisation insertions does not generate 
additional terms in the isovector kernel; a fact that very much simplifies the 
study of isovector systems, allowing one to employ a model gluon propagator 
and maintain Goldstone's theorem. 

These arguments are independent of the form of the model gluon propagator. 
The inclusion of 3— and 4— gluon vertices is a straightforward extension, which 
does not modify the conclusions. 

The isoscalar channel receives additional contributions from vacuum polari- 
sation diagrams. This may provide the basis for understanding the rj-r]^ mass 
splitting in the present framework[8]. 

^ In deriving this equation we have used the identity 

J dH5\lf-^f{k,p,l)^\ J dH5\l)p-kf{k,p,l) 
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Fig. 1. The replacement procedure of Eq. (10) illustrated to O(G^), which provides 
the kernel of the meson and diquark Bethe-Salpeter equations. The internal solid 
lines represent dressed quark and gluon propagators. 

The full 0{G^) kernel for the isovector-meson Bethe-Salpeter equation is illus- 
trated in Fig. 1. Using Eq. (8), the Bethe-Salpeter equation for isovector 
mesons is 



^M{p■,P) = -GJ^XMl^^ (11) 
- g 7^ ('5+ 'ju S+ 7^ xm + S+-i^XMl„S_ + xm Iv S- 7/. S_) 7^ 

where S± = S{p ± -P/2) and xm = S+ ^m{p', P) S-. Using Eq. (5) the Bethe- 
Salpeter equation for (3)c, isovector diquarks can be written 



^Up-^p)--\gx%i,-YqG'i, {s^i.s^^,x% (12) 

+5 S+ 7^ x% 7m S_ + x% Iv S- 7^ 5_) 7,. , 

with Fj-^s = r^gC* and x%3 = S+^d^^P' ^) '^h^re C = 7274 is the charge 
conjugation matrix. In these equations the underlined term is the "crossed- 
box" contribution of Fig. 1. 



Diquark bound- states. 

At 0{G) the only difference between Eqs. (11) and (12) is a 50% reduction 
in the coupling strength; i.e., the coupling is twice as strong in the meson 
equation. This follows from the algebra of the SU (3)c generators and entails [9] 
the existence of scalar and pseudovector diquark bound-states with m^^ > 
mL and ml^ > ml 



"qq "-"^ ■■"qq " ""qq 
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Of the terms in Eqs. (11) and (12) only the underhned one in each is repulsive 
in the 0^^, 1^^, 0^^, 1^^ channels. The numerical factors arise from the algebra of 
the SU (3)c generators and do not depend on the form of the gluon propagator. 
Relative to its companion, parenthesised contributions, this term is 5-times 
as large in the diquark equation, which provides for the possibility that it 
eliminates diquark bound-states. 

4. Results and Conclusions 

In order to study the bound-state spectrum one must solve Eq. (9). Our results 
are illustrated in Fig. 2. The O(G^) corrections become noticeable for <^ 
rf/l, which is important because it is the domain sampled in the Bethe- 
Salpeter equation. It is also important to note that one has dynamical chiral 
symmetry breaking for any G > 0; i.e., for rUq = 0: as{0) > 0, VG > 0. Further, 
(jy and as have no pole on the real-p^ axis. Therefore the quark propagator has 
no Lehmann representation and can be interpreted as describing a confined 
particle. 




Fig. 2. We plot crv{x) and crs{x) obtained as the solution of Eq. (9); 
(^viP^) = &v{x)/rj^, crs{p^) = &s{x)/r}. 

The model gluon propagator of Eq. (8) entails that the bound-state con- 
stituents have zero relative momentum; i.e., p = 0. In this case the most 
general form of the 0^^ Bethe-Salpeter amplitude is 



75 ■ 



(13) 
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Substituting Eq. (13) into Eq. (11), with p = 0, one obtains a 2 x 2 matrix 
eigen-value problem of the form H Q = Q; where the elements of H depend 
on P^, and 0"^ = {6i,62)- The cigcn- value problem is solved when that 
is found for which det{H{P'^) — /} = 0; the eigen- vector follows immediately. 
The same procedure is followed in each of the other channels. The vector- 
meson Bethe-Salpeter amplitude has the form 

(P) = ■ {P') + '-<^,^4pA' (P') , (14) 

with e^(P), A = 0, ±1, the polarisation vector: ■ P = 0; r^g*^ has the same 
form as the pseudoscalar meson amplitude in Eq. (13); and E]^'" the same 
form as the vector- meson amplitude in Eq. (14). 

Our chiral limit results are presented in Table 1. The 0{G):0{G) eigen- vectors 
are (^r,^^) = (0.83,0.55), (91,91) = (1,0). At 0{G^y.O{G^) they arc 
(91,91) = (0.87,0.49), (^{",^2") = (0-99, -0.12). 

Table 1 

Calculated meson and diquark masses, in GeV (?? = 1.06 GcV). The labels "0(G"'): 
0(G"^)" mean that the solution of the 0(G") quark Dyson- Schwinger equation was 
used in the 0(G"') Bethe-Salpeter equation. The non-zero quark mass was chosen 
to reproduce the experimental ratio rriT^/mp at 0(G^):0(G^). "Unbound" means 
that there is no solution of the associated homogeneous Bethe-Salpeter equation for 
real P^. 



rriq = 





0(G): 0(G) 


0(G): 0(G2) 


0(G2): 0(G2) 


m% 







0.30 









1.19 


Unbound 


Unbound 






0.750 


0.745 


0.823 


1+ 




1.30 


Unbound 


Unbound 


ruq = 


0.012 








^qg 




0.140 


0.328 


0.136 


0+ 

-rnqq 




1.21 


Unbound 


Unbound 


ml~ 




0.767 


0.760 


0.770 


1+ 




1.32 


Unbound 


Unbound 



The Goldstone theorem is manifest when the quark DSE and pseudoscalar me- 
son BSE are truncated consistently; i.e, at 0{G):0{G) and 0{G^):0{G^). This 
can be shown analytically. If the dressing is inconsistent; e.g., 0(G) :0(G^), 
the pseudoscalar is half as massive as the vector meson. 

The O(G^) corrections only provide for a small (10%) mass increase in the 
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vector- meson channel, as one would expect of a weak, net-repulsive interac- 
tion. It is weak because of the cancellation between the vertex-correction and 
crossed-box contributions, which is a necessary consequence of the preserva- 
tion of Goldstone's theorem. 

In the diquark channel, however, where the coefficient of the repulsive term is 
larger and the cancellation incomplete, the 0{G'^) corrections have the signif- 
icant effect of eliminating the bound-state pole on the real-P^ axis. 

We note that the pseudovector, ^2~, component of the pseudoscalar meson 
is a significant part of its Bethe-Salpeter amplitude whereas the tensor, $1 , 
component of the vector meson is small. This feature is also observed in the 
separable Ansatz studies of Ref. [10]. 

Our results for 7^ are presented in Table 1. The 0(G^):0(G^) meson 
eigen-vectors are: {9^', 9^') = (0.86,0.51) and (91', 9^') = (0.99,-0.13). The 
Goldstone boson character of the pseudoscalar is clear; i.e., at a consistent 
level of truncation its mass increases rapidly from zero as rUg is increased. In 
contrast, the vector mass shifts upward by < 1%. (The upward shift is ~ 10% 
for niq — O because the difference between the 0(G) and O(G^) solutions of the 
quark DSE is greater in this case.) A non-zero current-quark mass introduces 
little quantitative and no qualitative change in the diquark channels; i.e., they 
remain unbound. 

For ruq = 0.012 GeV we have plotted det{H{P'^) - 1} in Fig. 3. It illustrates 
the effect of the O(G^) repulsive term in the Bethe-Salpeter kernel, which 
shifts the zero in the meson channel very little but completely eliminates it in 
the diquark channel. The additional repulsive strength in the diquark channel 
is amplified by the growth in asix) for x < 0; i.e., the feature of confinement 
manifest in the quark propagator, and due to the infrared enhancement of the 
gluon propagator, plays a role in destabilising the diquarks. 

A motivation for our study is the observation that while the absence of a 
Lehmann representation for gluon and quark propagators ensures there are no 
gluons and quarks in the strong interaction spectrum, it neither entails nor pre- 
cludes the existence of bound-states, whether coloured or not. We have shown 
that the rainbow-ladder truncation is peculiar in the sense that it alone has 
no repulsive terms in the BSE kernel. In every truncation beyond this there 
are both attractive and repulsive terms. We saw that in a consistent trun- 
cation there was almost complete cancellation between the O(G^) attractive 
and repulsive terms in the meson channel. This is an indication of why the 
studies of meson spectroscopy and decays using a model gluon propagator in 
rainbow-ladder approximation have been successful. In the diquark channel 
the situation is quite different. The algebra of SU{3)c ensures that the coeffi- 
cient of the O(G^) repulsive term is larger and this, coupled with confinement 
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Fig. 3. det{H{P'^) — 1} plotted as a function of P^. This function vanishes at the 
square of the bound-state mass in the channel under consideration. 

as manifest in the form of the quark propagator, entails that the repulsive 
effect survives to ensure the absence of a stable diquark bound-state. This 
effect cannot reasonably be reproduced in rainbow-ladder approximation. It 
is plausible that these features persist at higher order and with the inclusion 
of explicit 3— and 4— gluon vertices. 
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